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Abstract. We have studied the dissipative dynamics of a solid-state qubit with an 
extra electron confined to either one of two coupled quantum dots. Previous theoretical 
work based on Bloch-type rate equations gave an unphysical uniform occupation 
probability of the electron in the quantum dots even for non-identical dots. We show 
that this is due to neglecting higher order interactions in the analysis. By including 
higher order terms, we obtain expected asymmetric occupation probabilities for non- 
identical dots. Our work demonstrates that the high order interaction terms can lead 
to important qualitative impacts on the operation of the qubit. 



PACS numbers: 73.2LLa, 73.23.-b 
1. Introduction 

Quantum measurement has attracted much attention because of its fundamental 
importance as well as its relevance to quantum computing. However, there are still 
vastly different points of view on quantum measurement. On one hand, the Zeno effect 
predicts that continuous measurement can freeze the state of a quantum system [1]. On 
the other hand, interaction with the outside environment such as a detector is expected 
to lead to the collapse of the state of the measured system. Therefore, reconciling these 
views by detailing the real mechanisms for selected model systems is very important. 

Many experiments on two-level quantum systems (see, e.g., Refs. [21 El Hj) have 
been performed to study this problem. Together with experimental investigations, some 
theoretical studies were also devoted to this important topic. Gurvitz et a/ [S, E] have 
theoretically studied a mesoscopic qubit-plus-detector system. The qubit consists of two 
coupled-quantum dots (CQD's) while the detector is a quantum point contact (QPC) 
interacting locally with one of the dots. The authors derived a set of Bloch-type rate 
equations to describe the measurement effects. Goan et al [7\ considered the same 
setup and developed a Lindblad quantum master equation using the quantum trajectory 
approach for dissipative processes. However, while reasonable result was obtained 
for a pair of identical quantum dots in both Refs. [5] and [7], these works predicted 
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unreasonable equal occupation probabilities in the dots when their energy levels are 
different. Subsequently, Li et al [8] developed a unified quantum master equation with 
the Markovain approximation and obtained distinct occupation probabilities for non- 
identical quantum dots as expected. 

In this work, we further study this qubit-plus-detector system following the original 
approach of Gurvitz et al [5]. We have found that the unphysical result of identical 
occupation probabilities in non-identical dots mentioned above is only due to neglecting 
higher-order terms in their derivation. We will show that higher-order terms are 
of fundamental importance. Retaining additional terms in our calculations leads to 
physically valid results consistent with those of Li et al [8j . Our analysis also illustrates 
the detailed mechanism of quantum measurement. 

2. Bloch-type rate equations 

The system we study is shown schematically in Fig. 1. It can be modeled by the 
Hamiltonian H = Hcqb + -f^QPC + -f^int, with 

HcQB=Eia\ai + i?24«2 + ^0(0102 + 4«i)) 
HQPc='^Eialai + ^^^0^0^ + '^Qir{a}ar + H.c), 

I r l^r 

Hint =^5Vtira\ai{a\ar + alai). (1) 

l,r 

Here, Hcqi) is the Hamiltonian of the CQD's. The Hubbard terms are neglected because 
we consider only one extra electron in the CQD's. The parameter Qq denotes the hopping 
amplitude of the extra electron between the two single-dot states (see Fig.l). Hqpq is 
the Hamiltonian of the QPC, and ifjnt is the interaction Hamiltonian between the CQD's 
and the QPC. Electrostatic effects are included in Hi^t. an extra electron in the left dot 
will lead to an effective variation 6Qir in the coupling Qir between the states Ei and Er 
in the two reservoirs of the QPC. The resulting coupling then becomes Q'l^ = Qir + 6Qir. 

In the occupation representation, the many-body wave function of the entire system 
(the CQD's plus the QPC) can be written as 

\^{t)) = [bi{t)a\ + b2it)al+J2biir{t)a\alai + J2biri{t)a\alar 

l,r l,r 

+ ^b2ir{t)alalai + '^b2ri(t)alajar + ^ biii>rr'(t)a\alal,aiau 

l,r l,r i<l',r<r' 

+ X! birrni'{t)ciio\a\,arar' + ^ b2Wrr'{t)<Aalal,aiaii 

Z<l',r<r' l<l',r<r' 

+ b2rr'ii'(t)alalalarar' ^ ]|0), (2) 

l<l',r<r' 

where bj(t), j = 1,2, 1/r, 2/r, • • -, are the time- dependent probability amplitudes of 
finding the system in the corresponding states. The vacuum state |0) corresponds to 
the state of the entire system with no extra electron in the CQD's and before the transfer 
of any electron between the reservoirs in the QPC. Then, bur{t) for instance denotes the 
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Figure 1. Two coupled quantum dots (CQD's) connected with a quantum point 
contact (QPC). Initially, all states up to the Fermi energy level {^l or fin) in each 
reservoir are filled by electrons. Ei and E2 denote the energy levels of an extra electron 
in the left and right quantum dots, respectively. The location of this extra electron 
will affect the current through the QPC. The quantum state of the CQD's can hence 
be read out by a current measurement. 



probability amplitude that the extra electron is in the left dot and an electron has passed 
through the barrier in the QPC. The initial conditions are &i(0) = 1 and bj{0) = 
for j = 2, llr, 2lr, ■ ■ ■. Substituting the many-body wave function [Eq. ([2])] and the 
Hamiltonian [Eq. ([T])] into the time-dependent Schrodinger equation ili'it)) = H\ip{t)) 
gives an infinite set of linear equations for the probability amplitudes hj (t) . Performing a 
Laplace transform to each equation, the transformed amplitude bj{E) = Jf^°° e^^%j{t)dt 
follows: 



{E-E^)h^{E) 


-QobiiE)- 


-Y^QlburiE) ~Y.^[,^Mr'v{E) = 2, 

l,r l',r' 

- ^^lrb2lr{E) — ^Ql'r'b2r'l'{E) = 0, 
l,r l',r' 






(3 a) 


{E-E2)hiE) 


-QobiiE)- 






m 


(E + El — El - 


■Er)biir{E) 


— Q'ij.bi{E) — Qob2lr{E) — Q',,^,biu'rr 

I'r' 


(E) 


= 0, 


(3c) 


{E + E1-E2- 


-E,)b2lr{E) 


— ^lrb2{E) — Qobiir{E) — ^l'r'b2U'rr 

I'r' 


(E) 


= 0, 


{3d) 


(E + Er-Ei- 


-Ei)birl{E) 


— Q'i.^.bi{E) — Qob2rl{E) — Q!i,,,.,birr'll 

I'r' 


(E) 


= 0, 


(3e) 


(E + Er — E2 - 


-Ei)b2rl{E) 


— Qirb2{E) — QQbirl{E) — y^ ^l'r'b2rr'll 


(E) 


= 0, 


m 



I'r' 



To solve Eqs. (3), Gurvitz et al. [5l|6] introduced an approximation by keeping terms up 
to the order of Vt^ to simplify the corresponding sum terms [see Eqs. (A5) in Ref. 5]. After 
tracing out the degrees of freedom of QPC, a coupled set of Bloch-type rate equations are 
derived. As shown in Fig. 3 of Ref. 5(b), however, the resulting occupation probabilities 
for states \Ei) and \E2) are 0.5 when t — > 00 for both symmetric and asymmetric CQD's. 
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This is correct for the symmetric case with identical quantum dots, but is unphysical 
in the latter case. Here, we show that this unphysical feature can be rectified by taking 
terms of the order of VL^VLq/Vci or even higher into account. Also, following Ref. 5, we 
consider the high voltage regime with eVd 3> where Vd = Hl — f^R is the apphed 
gate voltage, and p is the density of states in the reservoirs. In this high voltage regime, 
when the temperature is low enough (here it is chosen zero for simplicity, as in Ref. 5), 
the electron has an extremely low probability to pass through the QPC from the right 
reservoir (with a lower pn) to the left one (with a higher /i^). Thus, the terms describing 
the back processes, e.g., J2i<v,r<T:' ^irr'wit), can be neglected. 

Taking the sum term Y^i ,^Q'i^biir{E) in Eq. (!3aj) for example, we obtain that 

W 7 A 

E QlhUE) ^-—{Vd + E- E,)h{E) + ^h{E), (4) 

where D' = 2'KpiprQ''^Vd, A = 2'KpiprQoQQ'Vd, and the corresponding exact calculation 
are shown in the Appendix. With the same approximate treatment, the resulting 
equations for bj{E) are 

[E-E^ + '^{Vd + E- E,)ME) - (no + ^ME) = t, (5a) 
ZVd ZVd 

iD iA 
[E-E2 + —{Vd + E- E2)ME) - (1^0 + t^ttME) = 0, (56) 

iD' 7 A 

[E + Eiri + —{Vd + E + Ei,,)]hur{E) - n'h{E) - {Qq + —)b2lr{E) = 0, (5c) 
-^Vd ^vd 

[E + Eir2 + ^(l^d + E + Eir2)]h2lr{E) - ^h2{E) - fl^o + 7^) hlr{E) = 0, {hd) 

ZVd \ zVdJ 



where D = 27rpiprfl'^Vd and Eirm = Ei — Em — Er {m = 1,2). Now, we introduce the 
notation 

For instance, crj^'' = J2biir{t)bli.^{t) denotes the occupation probability for the extra 

Ir 

electron staying in the left quantum dot and an electron passing through the QPC. 

The equations for the amplitudes bj{E), j = 1, 2, 1/r, 2/r, ■ ■ -, can be converted to 
a new set of equations for a^j'' using the inverse Laplace transform: 



l-r- 



We now multiply Eq. fl5cj) by h\i^.{E'). The resulting equation is then subtracted by its 
complex conjugate after exchanging E with E' . We obtain 
dEdE' 



4^2 Y.{iE' - E - iD')bur{E)hl,AE') 

Ir 

'-^[{E' + + {E + E,,{)]hur{E)hl,^{E') 



+:j^[b2UE)blUE') + bur{E)b;UE')] - no[buriE)b;UE') - b,lr{E)b\,,{E')] 

-n'[bl{E')bur{E) - 6i(E)6*,,(E')]}e^(^'-^)* = 0. (6) 

Using Eqs. 4, fl5aj) and fl5cj) . we can obtain from Eq. ([6]), by neglecting terms of order 
0{n^/V^) and higher, that 
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Similarly, Eq. ( l5dj) gives 
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4ft + 
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2V,, 
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Ol2 021 J 



cr. 



(1) 

21 



Ol2 



4?] 



(7) 



(8) 



To calculate the off-diagonal element au, we subtract the product of Eq. (l5q) and 
b^riE') by the product of bur^E) and the complex conjugate of Eq. fISdj) after exchanging 
E ^ E'. We obtain 



A 



A 
2^ 



^2^ =^^^12^ " i^^' + + V^(Ti°2^ + ^^]ok^}^ - ^22 



(l)l 



A 

2V^ 



^' (0) , ^ (0) 



41^ + 



(1) 

22 



+ 



A 



(1) 

11 



(1)1 
22 J) 



(9) 



where e = E2 — Ei is the energy-level difference of the two quantum dots. Similar 



procedures can be used for 4^^ > 2) and the resulting equations are 

^(")_ n'[^(") ^("-1)1 I VO rn-(") -U ^ f^'^ Fn-^") 4- n- 

(Tn ---^^ [0-11 - (111 J+ZiZo[(Ji2 - 0-21 J + (^i - [(J12 +(^21 " 0" 



(n-1) 
12 



(n-l)n 
21 J 



<y22 ■ 



Ol2 ■ 



An) 



2Vd 



An) 
^2 



in) (n-1) 
^21 - ^12 



4rft 



- \{D' + D)42^ + VDD^air'^ + zno[4-^ - aS>] + ^[4?^ + 4?] 



A 



n' (n-1) 

TT^ii 



f^^("-i) 

n^"22 



+ ^f4r'^+4rft 



f2^(n) I O' (n) 

syOii + o" 022 



A n (n) n' (n) 
2V^ril T-'^22 J - 2V^ [sy^ll + "0 ^22 J • U^J 

Summing over n, we obtain the following Bloch-type rate equations for the density- 
matrix elements: 



&ii=iQo{ai2 - 0-21), 

(T22=ii7o(o"21 - 0-12), 

ai2=i€ai2 - -^'^i^ + ^^^0(0-11 



0"22) 



X. 



0"11 + 0-22j 



with 



(12) 



Here cTj,- = J2n 4^^ is the reduced density-matrix elements of the CQD system after 



tracing out the variables of the QPC The decoherence rate characterizes the 
exponential damping of the off-diagonal density-matrix element. These Bloch-type rate 
equations are the same as those from a theory based on master equations with the 
Markovain approximation [cf., Eq. (17) in Ref. 8]. When ignoring the effects of the 
higher-order terms x = 0, these equations are reduced to the rate equations derived 
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Figure 2. (Color online) Occupation probability of the extra electron in each quantum 
dot as a function of time for = Hq and (a,) e — E2 — Ei = 0, x = 0; (b) e = 0, 
X = O.mo; (c) £ = r!o, X = 0; and (d) e = ^o, x = O-l^o- 



in Ref. 5(b), which gives the unphysical result of identical occupation probabilities for 
asymmetrical CQD's. This implies that the higher-order terms play an essential role 
in obtaining more accurate and reasonable results and help us better understand the 
effects of the quantum measurement on the considered system. 

Here we take into account the terms of the order fl'^flQ/Vd, while these terms were 
ignored in Ref. 5. To obtain more accurate results, one can consider further higher-order 
terms, which have the order of Q^/V^ and Q^Qo/V^. In the high voltage regime, these 
higher-order terms do not affect the results significantly (see the Appendix). 

3. Discussion and conclusion 

Figure 2 shows the time dependence of the occupation probability of the extra electron 
in each quantum dot. When x = 0, the occupation probability in each dot always 
decays to 0.5 at t 00 [see Figs. (2a) and (2c)]. These results are identical to those in 
Ref. 5(b). However, after including higher order terms, x becomes nonzero and is not 
negligible. This gives rise to an additional term ^x{(^n + o"22) in the last equation of 
(11). From Fig. (2b), one can see that the occupation probabilities still decay to 0.5 as 
t 00 for symmetrical CQD's, but decay more slowly than in the previous case with 
vanishing x- Moreover, for asymmetrical CQD's, the occupation probability in each dot 
decays to a different value at t — >• 00 [see Fig. 2(d)]. This is reasonable because the two 
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quantum dots are characterized by different parameters. Below we further explain these 
phenomena analytically. 

Let us look at the stationary solution of Eq. (11) in the limit t — >■ oo so that 
&ij(t — > oo) = 0. Using aij{t oo) = Uij + ivij, where Uij and are real numbers, 
Eq. (11) can be rewritten as 

0= - 2noVu 
0=21]ot^i2 

0=ie{ui2 + ivu) - y(^i2 + if 12)+ «^o(crii - (^22) - ^ • (^3) 
From these equations, one obtains 

(^U - = TZ-^ = TT- (14) 

As expected, an = (J22 for Ei = E2 (i.e., e = 0), while an 7^ a22 for asymmetrical 
CQD's with El ^ E2 (i.e., e ^ 0). 

In conclusion, we have presented a quantitative description of the dissipative 
dynamics of a CQD qubit connected to a detector in the form of a QPC. Bloch-type rate 
equations for the reduced density matrix elements of the qubit are derived. In particular, 
the long-time probability distribution of the state of the qubit is found to depend on the 
energy levels of the single dot electron states. This corrects an earlier calculation which 
has predicted a distribution independent of the dot properties. The improvement results 
from taking into account higher order interaction terms in our analysis and it can be 
extended to other quantum systems. In the present work, the Bloch-type rate equations 
are derived for a CQD system working at zero temperature and applied in the strong 
bias voltage regime. However, a finite temperature will affect the asymptotic population 
of the qubit and a popoulation inversion can even be achieved in some cases [9] . 

4. Appendix: Derivation of equations for probability amplitudes 

For the QPC, we assume that the hopping amplitude Qir{Ei, E^) between the left and 
right reservoirs depends only weakly on the energy levels Ei and Er. The energy 
dependence is therefore neglected, i.e. Qir{Ei,Er) = and Q'ir{Ei, E^) = fl'. As 
high-order terms J2i'r' ^'biii'rr'{E) and J2i'r' ^b2ii'rr'{E) are neglected, Eqs. (3c) and (3d) 
give 

{E + E1-E2- Er)n'%i{E) + n'nonb2{E) 

b,ir[t.) ^ ^ _ _ ^^^(^ + E,-E2- E,) - Ql 

Because the energy levels in each electron reservoir of the QPC are dense, we can 
replace each sum over / and r in Eqs. (3a) and (3b) by an integral. For instance, 
J2ir I Pi{Ei)pr{Er)dEidEr, wherc pi^r is the density of states in the left and right 
reservoirs, respectively. This integral can be split into two parts: the principal and 
singular value parts. Here the bandwidth of the QPC reservoirs are much larger than 
Vd and the principal part is thus negligibly small pjO]. Actually, the principal part 
merely renormalizes the energy levels and the singular value part plays the dominant 
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role. When Eq. (A.l) is substituted into J2ir^'biiriE), we thus obtain two terms. The 
first term is given by 

{E + El - E2 - Er)n'%,{E) 
j^{E + Ei-E,- Er){E + Ei~E2- E^) - fig 

/^^^^ /■+-^^ Pi{Ei)pr{Er)n'\E + El - E2- Er) 



dEi / dEr- 



'm {E + El — Er — E[){E + Ei — Ej. — E'2, 
hi{E) / dEi / dEr{-mpipr^ ) x 



E + El — E^ — El, 



{E + Ei-Er- E'^y + r/2 
E + El — Er — E[ 



2 5{E + El - Er - E[){E + El - E2 - Er] 
5{E + Ei-Er- E'2){E + Ei-E2-Er 



{E + Ei- Er - E[f + 7]^ J 

-^^fei(^) [\,{^ + \)iy, + E- E[)e{V, + E-E[) 

+^(^ 'l)iV, + E- E'2)e{V, + E-E'2)\, (A.2) 



where Va = pi - I^b. D' = 2npiprn'Wd, e = E2 - Ei, e' = ^ + EI2 = 
+ E2) =F and 9{x) is the Heaviside step function. In the high voltage hmit 
with eVd ^ fi^p, the step function in Eq. (A.2) becomes one. Thus, the L.H.S. of 
Eq. (A.2) can be approximated by —i{D'/2Vd){Vd + E — Ei)bi{E). When higher-order 
terms, i.e., the terms of the order 0{Q^/V^), are inculded, Eq. (A.2) is approximated 
by -i{D' /2Vd){Vd[l - ^{D'/Vdy] +E- Ei}bi{E) in the high voltage limit. Obviously, 
the correction by {D' /Vd)^ is very small for a large voltage Vd- 
The second term is 

\f{E + Ei~E^- Er){E + E1-E2- Er) - 

/^^^^ /•+-^^ Pl{Ei)pr{Er)n'QQo 



dEi / dE, 

-00 Jria 



IJ.R (E + El — E'l — Er){E + El — E2 — Er) 

r^iL r+co ( \ 



j-liL r+co I 

h2{E) / dEi / dEri-i^piPr^'nn^) — — -6{E + Ei-E'^-Er 

-6iE + Ei~E[-Er)] 



E + Ei- E[ 
1 



E + El — E'2 — Er 
—iA 1 

= ^^2(i?);^rr^ liya + E~ E'2)9{Vd + E-E'2)-{Vd + E- E[)e{Vd + E- E[)] , 

(A.3) 

where A = 2tt piPr^'^lo^lVd- In the high voltage limit, the two step functions in Eq. (A.3) 
become one and the L.H.S. of Eq. (A.3) is finally approximated by i{A/2Vd)b2{E). Thus, 
the sum of these two terms gives Eq. (jl]) in the high voltage limit. Substituting Eq. (jlj) 
into Eq. (3a), one obtains Eq. ( l5a|) . Also, one can derive Eqs. ( l5B -(l5dl) using similar 
procedures. When higher-order terms, i.e., the terms of the order 0{Q^Qo/Vj), are 
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considered, Eq. (A. 3) is approximated by 2(A/2V^/)(1 + j)b2{E) in the high voltage 
regime, where the correction 7 = [{D + D')yDD' — DD']/AV^ is very small for a large 
voltage Vd- 
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